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ABSTRACT 

Continuum emission at 220 and 355 GHz from the prototype Class source 
VLA 1623 has been detected using the JCMT-CSO interferometer. Gaussian 
fits to the data place an upper limit of 70 AU on the half-width half-maximum 
radius of the emission, which implies an upper limit of ~ 175 AU for the cutoff 
radius of the circumstellar disk in the system. In the context of existing collapse 
models, this disk could be magnetically supported on the largest scales and 
have an age of ~ 6 x 10^ yr, consistent with previous suggestions that Class 
sources are quite young. The innermost region of the disk within ~ 6 AU is 
likely to be in centrifugal support, which is likely large enough to provide a 
drive for the outflow according to current theoretical models. Alternatively, if 
175 AU corresponds to the centrifugal radius of the disk, the age of the system 
is ~ 2 X 10^ yr, closer to age estimates for Class I sources. 

Subject headings: accretion: accretion disks - circumstellar matter - stars: 
formation - stars: low-mass, brown dwarfs - stars: pre-main-sequence 
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1. Introduction 

Submillimeter observations of deeply embedded star-forming regions have identified 
a small number of strong sources that have virtually no emission at wavelengths below 
10 lira and whose spectral energy distributions are characterized by single blackbodies 
at T ~ 15 — 30 K. These cold "Class 0" sources are good candidates for protostars 
whose circumstellar envelope masses exceed that of their central protostars (Andre, 
Ward- Thompson, & Barsony 1993, Andre & Montmerle 1994, Andre 1995), in which case 
their ages may be as low as 2 x 10^ yrs (Barsony 1994). Sources in this class are also defined 
by the presence of energetic bipolar outflows. 

There is now good evidence that a large fraction of young stellar objects have dense, 
dusty accretion disks (eg. Beckwith et al. 1990, Lada 1991). This is in accord with standard 
collapse models of star formation (eg. Terebey, Shu, & Cassen 1984), which predict that 
systems older than 10^ years have centrifugally supported disks extending to 100 AU. 
Disks around the sources L1551 IRS5 and HL Tau have been resolved with semi-major 
axes of 80 and 60 AU, respectively (Lay et al. 1994, Mundy et al. 1996). Even larger 
magnetically supported "pseudo" disks extending out to 1000 AU are predicted for some 
magneto- hydrodynamic collapse models (Galli & Shu 1993, Fiedler & Mouschovias 1993). 

Most current models of outflow posit the existence of disks as an essential element of 
the engine (eg. Pudritz, Pelletier, & Gomez de Castro 1991, Konigl & Ruden 1993, Shu 
et al. 1993). If the Class sources are young, then the Terebey et al. (1984) models 
suggest that a disk with radius of 100 AU at 2 x 10^ yrs has a centrifugal radius of only 
Tc oc ~ 0.8 AU at 2 X 10"^ yrs. Similarly a larger magnetically supported outer region of 
rs oc f/^ ~ 20 AU (GaUi & Shu 1993) should be present at these early times if 1000 AU 
structures exist at 2 x 10^ yrs. Such small, centrifugally supported regions would still be 
able to drive outflows, while the larger magnetic structure could play an important role 
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in star formation itself. This Letter describes the results of our search for an accretion 
disk in the prototype Class source, VLA 1623, using the JCMT-CSO submillimeter 
interferometer. 

2. Observations and Analysis 

Observations of VLA 1623 were obtained with the JCMT-CSO interferometer on 7 
April 1995 at a frequency of 356 GHz with 500 MHz bandwidth and 9 April 1995 at 353 
GHz with 1 GHz bandwidth. The projected length of the baseline ranged from 140 to 
190 kA and observations were obtained over a period of two to three hours each night. 
The quasar NRAO 530 was used to monitor the gain of the instrument and the absolute 
uncertainty in the flux calibration is 20%. The average gain of the instrument was 70 Jy 
K^^ and the system temperature was 800 K. The data were obtained with 10 s integrations 
and vector-averaged over 100 s intervals. Since the phase response of the interferometer 
over longer periods of time is not yet known, only the amplitude of the visibility function 
was used in the analysis. Observations were also obtained at a frequency of 220 GHz with 
1 GHz bandwidth on 10 April 1995 over a period of 7 hours. The projected length of the 
baseline ranged from 60 to 120 kA and the system temperature was 400 K. 

The visibility flux of VLA 1623 shows no dependence on the projected baseline length 
(Figure ^,b) and thus the source is unresolved. We measured both the flux and an upper 
limit to the full-width half-maximum radius of the source using the techniques outlined in 
Lay et al. (1994). Briefly, we fit a two-dimensional gaussian to the data using the Rice 
distribution to describe the noise. The position angle was fixed at 35° to orient the disk 
perpendicular to the large-scale outflow (Andre et al. 1990). Fifty simulated data sets were 
then generated by adding artificial noise to the best fitting model. The best fit parameters 
(flux and major-axis radius) derived from these simulated data sets are shown in Figure [^d. 
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from which 90% confidence hmits can be determined. At 355 GHz, the fiux is 0.5lo;i Jy and 
the major-axis radius is -R < 70 AU (for a distance to p Oph of 160 pc). Fitting the data 
with a circularly symmetric gaussian produces a somewhat smaller upper limit to the radius 
of 45 AU. A more physically valid model would be a power-law disk (see §3). The large 
uncertainty in the fiux at 355 GHz is due to the rather limited baseline coverage obtained 
for this low declination source. As a result, the model fits are unable to distinguish between 
a weaker compact source and a stronger, more extended source for which some of the fiux 
is resolved out even on the shortest baselines at 355 GHz (Figure |ID. At 220 GHz, the fiux 
is 0.12t[J:[]^ Jy and the radius is R < 150 AU. 

We calculated the expected contribution to the visibility fiux from the extended 
envelope as in Lay et al. (1994). We modeled the envelope of VLA 1623 as an optically thin 
sphere with T cx r~°-^ and p oc r^^-^. The fiux scale was calibrated using single dish fiuxes 
from Andre et al. (1993). We also assumed that the envelope stops where the disk begins 
and so adopted an inner cutoff radius for the envelope of 70 AU. Under these conditions, 
the envelope for the 70 AU model contributes negligible flux at 355 GHz and only 0.03 Jy 
at 220 GHz. Thus, we attribute most of the observed flux to the presence of an unresolved 
disk in VLA 1623. 

3. Fitting Accretion Disk Models 

Our observations provide useful constraints on plausible accretion disk models for the 
origin of the submillimeter emission. We adopt the approach of Beckwith et al. (1990) 
and assume a geometrically thin disk with power law scalings with radius of the disk 
temperature and surface density, T(r) = Toirjrf^'^ and S(r) = Y,o{r j r ^ where Tq and 
So are the dust temperature and mass surface density at the inner radius of the dust 
distribution. The optical depth is r = Hn^/cosd, where 6 is the inclination of the disk axis 
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to the line of sight and = >^oiyl'^o)^ is the specific dust opacity. The fiux emitted by the 
disk is then = {27rcose/D^) l^f B^{T){1 - e'^ydr where B^(T) is the Planck function, 
D is the distance to the source, and Rd is the outer disk cutoff radius. Power-law fits to 
disks around HL Tau and L1551 IRS5 produce cutoff radii up to a factor of 2.5 larger than 
radii from gaussian fits (Lay et al. 1994, Mundy et al. 1996). However, since VLA 1623 is 
unresolved, there is no structure in the visibility data to constrain a power law and so such 
a fit was not performed on this data set. 

While q = 0.75 for purely viscous accretion disk models, more general calculations 
which incorporate the effects of reprocessing and disk fiaring typically find 0.5 < g < 0.75. 
Noting that various model fits have 0.5 < p < 1.5 (eg. Thamm, Steinacker, & Henning, 
1994), we adopt ranges of 0.25 < q < 0.75 and -0.25 <2 - p - q < 1.25. The maximum 
dust temperature must be less than its sublimation temperature, Tg^max ^ 2000 K. The 
inner edge of the disk cannot be inside the stellar photosphere, which for low mass stars 
gives To > -R* ~ 2 = 0.01 AU. We take the inclination angle of the disk to be that of the 
outflow from VLA 1623, 9 = 80° ± 5° or almost edge-on (Andre et al. 1990). The spectrum 
oc z/" determined from our measured fluxes at 355 and 220 GHz has a ~ 3.01212) where 
optically thin models have a — 2 + P and optically thick models have a — 2. Recent 
calculations show that /3~1.5±0.5at wavelengths from 650 //m to 2.7 mm for submicron 
to millimeter sized particles and a wide variety of densities and temperatures (Pollack et al. 
1994). Finally, the parameter space characterizing the emission is also constrained by our 
observed flux, our upper limit to the cutoff radius of ~ 2.5 x 70 = 175 AU, and by the ~ 1 
Lq bolometric luminosity of the entire source (Andre et al. 1993). 

It is obvious that this accretion disk model covers a large parameter space that is 
relatively poorly constrained. The issue of the uniqueness of model fits has been addressed 
by Thamm et al. (1994) and Bouvier & Bertout (1992), who find "global ambiguity" 
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in acceptable fits to models of T Tauri disks based on spectral energy distributions that 
incorporate far more infrared data than is available for VLA1623. However, for the limiting 
cases of completely optically thick or thin truncated disks in the Ray leigh- Jeans limit, the 
cutoff radius can be related to the observed flux Fi, by 

Ro/r^ = [{^Ljxcose) + 1]^/^ (1) 

where Ly = D'^Fy/2T[{2kv'^ / (?) can be computed from the observations. For the optically 
thick case, the parameters are 7 = 2 — g and x = ToT^, while for the optically thin case, 
7 = 2— p — g^^O and x = ToTotI, where Tq = HyT^o/cosO. 

Equation (1) shows that we can simplify our analysis of the model by reducing the six 
free parameters {Ky, q, p, To, Tq, Sq) to essentially the two parameters 7 and x (with Rd 
measured in units of Tq). Then a disk model characterized by a given index 7 produces a 
unique curve in the Rd/^o ~x plane. These curves are plotted in Figure ^for VLA 1623, for 
which = 0.5 Jy gives L355 = 2.25 x 10^ K AU^. In addition, our observed flux and upper 
limit to the cutoff radius restrict possible disk models to lie within the region demarcated 
by the three dashed lines in Figure |[ Boundary line (A) is the 7 = 2 model (i.e. g = or 
p + g = for optically thick or thin cases), which demarcates a model disk with uniform 
surface brightness and therefore the smallest cutoff radius Ro for the emitted flux. Line 
(B) is our upper limit to the cutoff radius Rd < 175 AU coupled with the condition that 
To > 0.01 AU, while line (C) is the upper limit to Rd/tq for a given x. Line (C) arises from 
the fact that To must be less than To^max, so that Tq must be greater than {x/To^max)^^'^ ■, 
and thus Ro < 175 AU implies RdI^o < 175(To „iaa./x)^/^. Figure ^ shows that while 
optically thick models with values of g ^ 0.25 — 0.5 fill the largest portion of the permitted 
parameter space, models with g up to 0.75 are also acceptable. Optically thin models fill a 
much smaller portion of parameter space, with low values of 2 — p — q being less likely. In 
general, the optically thin solutions have relatively large values for the inner radius, which 
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suggests that they may also have a small optically thick region that contributes relatively 
little to the total flux of the source. 

The observed bolometric luminosity of the envelope+disk system also provides a 
constraint on the cutoff radius. For q ^ 0.5, the luminosity of an optically thick disk is given 
by = A7iaTXi{RD/rof-^'^ - l)/(2 - 4g), while for q = 0.5, U^i = AnaTyMRD/r^). 
For q < 0.5, this expression can be combined with equation (1) to give 

R . . TT (2 - Q?^' ( Lbol A -1/6 ^ 1.355/2.25 X 10^ K AU^ ^2/3 
i?^ = 47.7AU^2-4g)i/6iTL;J i cose/cos80 ) 

In addition, equation (1) can be rewritten as ToV^ oc R'j^'^, which allows the straightforward 
calculation of Tor^ for a given Rd obtained from equation (2). For g = and L^oi < 1 Lq, 
the cutoff radius is Rd > 68 AU and the inner disk temperature is To < 57 K, while for 
q = 0.25 and Vo = 0.01 — 10 AU, the radius is >69 AU and the temperature is To < 430 — 80 
K. For g > 0.5 the expression for R^ is different from equation (2) and depends on the value 
of To- For q = 0.5 and = 0.01 — 1 AU, the cutoff radius is >90-85 AU, while for q = 0.75 
the radius is >230-120 AU. A smaller but warmer optically thin disk could probably also 
produce the observed 355 GHz flux and bolometric luminosity, but the numerical models 
required to derive the exact cutoff radius in this case are beyond the scope of this paper. 

Thus both the spectral index and the models summarized in Figure |^ are consistent 
with either optically thin or optically thick disks with a cutoff radius not exceeding 175 
AU. We obtain a lower limit to the disk mass from the measured 220 GHz flux by taking 
the optically thin limit, 

M > {D^F,/k,B,{To))[2-p- q/2-p]{Rn/roy (3) 

Assuming a cutoff radius Rd = 175 AU, an inner radius = 0.01 AU, an opacity k^, = 0.01 
cm^ gm"^ at 220 GHz {f3 = 1.5), a maximum temperature of To^max = 430 K, and q ~ 0.25 
gives a minimum disk mass M > 3 X IQ-^ Mq; the true disk mass could be much higher. 
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4. Discussion 

The constraints on the size of the central source have imphcations for theoretical 
models of the evolution of sources such as VLA 1623. We adopt an effective sound speed 
a = 0.27 km s~^ for this innermost region of VLA 1623 (Andre et al. 1993) and a magnetic 
field strength B — 30 /iG appropriate for low mass molecular cloud cores (Myers & 
Goodman 1988). In the context of the GaUi & Shu (1993) model, the magnetic support 
radius of the disk is given by 

TB = 0.12(6-2570)'/^/=^ = 2.75(5/30 /xG)^/=' (a/0.27 km 8'^^^%^^ AU (4) 

where ^4 is in units of 10^ yrs and G is the gravitational constant. If we identify the upper 
hmit to the cutoff radius of 175 AU with the radius of the pseudo-disk, rs, then this 
model predicts an age for the system of t < 5.9 x 10''(-B/30 /xG)~'^/^(a/0.27 km s"^)^/'' yr. 
This age agrees with general estimates for Class ages (Barsony 1994), but is somewhat 
longer than initial estimates for the age of VLA 1623 (Andre et al. 1993). In a 
collapsing singular isothermal sphere (SIS) model, the expansion wave has reached a radius 
T^exp — dtexp < 3400(tesp/5.9 X 10"^ yr) AU, which is comparable to the size of the 2000 
AU envelope seen in the Andre et al. (1993) map. In the context of the SIS model, this 
age suggests that the envelope and the collapsed protostellar core+disk will each contain 
comparable amounts of material. 

Alternatively, if we identify 175 AU with the centrifugal balance radius, 
rc = O.ObSaQlt^ = 0.030(f)o/10-^3g-i)2t| AU, then the age of the system is < 1.8 x 10^ yr, 
similar to age estimates for Class I sources (Wilking, Lada, & Young 1989). In this picture, 
the magnetic pseudo-disk, if present, would extend to 2300 AU, and must have a very 
low column density in order to be undetectable with lower resolution measurements (e.g. 
Andre et al. 1993). This age is five times longer than the time required for the expansion 
wave to reach the edge of the circumstellar envelope. Over this time, a compressional wave 
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propagating inwards from the envelope edge would have accelerated the collapse of the 
remaining envelope material onto the stellar core. 

Returning to the assumption that most of the compact core is magnetically supported, 
the centrifugal balance radius in the pseudo-disk then occurs at Vc < 6.2(f2o/10~^^s~^)^ 
AU. This region is large enough either to drive off a disk wind or to interact with a possible 
stellar magnetosphere, and thus theoretical models developed to explain outflows in Class 
I sources are likely to be equally vahd for Class sources. The total collapsed mass of the 
star+disk system in the SIS model is M^ou ~ [1 + {'K/A)Y^2a^rexp/G = 0.29 M©. Taking 
a stellar radius — 2.5 (vahd for models with < 0.3 Mq, Stabler 1988) and 
a bolometric luminosity L* ~ 1 (Andre et al. 1993), we find that an accretion rate 
through the disk of Ma = L^R*/GM^ = 0.27 x 10"^ M© yr"^ can provide the observed 
luminosity of VLA 1623. This accretion rate is significantly smaller than the rate inferred 
from the collapse of an SIS model, Mj = 0.975a^/G = 4.5 x 10~^ Mq yr~^ A discrepancy 
of a factor of 5 between the disk and collapse accretion rates in SIS models has been noted 
by Kenyon, Calvet, & Hartmann (1993) for more evolved T Tauri systems. A possible 
reason for the inconsistency in the mass accretion rates is that SIS models have constant 
infall rates. More general models of molecular cloud cores (eg. McLaughhn & Pudritz 
1996, 1997, Foster & Chevalier 1993) do not collapse at a constant accretion rate, which 
may alleviate this problem. 

Finally, we note that the outflow in VLA 1623 has reached much larger scales (~ 20000 
AU, Andre et al. 1990) than the expansion wave predicted by the SIS model. This 
situation is possible because outflows are likely to be super- Alfvenic. The effect that such 
an outflow has upon collapse has yet to be investigated. 



5. Conclusions 



High resolution submillimeter continuum observations have detected a compact source 
in VLA 1623, which is most readily interpreted as a disk. Gaussian fitting gives an upper 
limit of 70 AU for the half-width half-maximum radius, which implies an upper limit of 
~ 175 AU for the cutoff radius of the circumstellar disk in the system. In the context of 
existing collapse models, such a structure could be magnetically supported on the largest 
scales and have an age ~ 6 x 10^ yr. The innermost region of the disk within ~ 6 AU is 
likely to be in centrifugal support, which could provide a drive for the outflow according 
to current theoretical models. Alternatively, if the entire observed disk is in centrifugal 
support, then the age of the system is ~ 2 x 10^ yr, similar to age estimates of Class I 
sources. In this case, we may be detecting a Class I source viewed nearly edge-on. 
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behalf of the Particle Physics and Astronomy Research Council of the United Kingdom, 
the Netherlands Organization for Scientific Research, and the National Research Council of 
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Fig. 1.— (a) The visibility flux of VLA 1623 at 220 GHz detected with the JCMT-CSO 
interferometer as a function of the projected baseline length. The flux shows no dependence 
on projected baseline length and thus the source is unresolved, (b) The visibility flux at 355 
GHz as a function of basehne length, (c) The visibility flux at 355 GHz as a function of the 
hour angle of the observations, (d) The source flux at 355 GHz deduced from 50 Monte Carlo 
simulations as a function of the deduced upper limit to the source size. With the available 
data it is impossible to distinguish between a weak compact source and a more extended, 
somewhat stronger source. The dashed line indicates the 90% confidence level upper hmit 
to the radius of the disk. 

Fig. 2. — (a) Cutoff radius Rd of an optically thick disk in units of the inner disk radius 
versus the product T^r^, where is the temperature of the disk at radius To- Solid lines 
are model fits for different values of g, the exponent in the temperature gradient. Viable 
accretion disk models lie within the region demarcated by the dashed line segments, A,B,C, 
described in the text. The angle between the disk axis and the hue of sight is ^ = 80°. (b) 
As in (a), but for optically thin emission. In this case Rd depends on the product ToTgrl. 
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